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FOREWORD 


The  nuterlal  preiemed  In  thii  Teport  appeared  earlier  as  a  portion  of  USL  Report 
No.  513,  Colottus  11  Summary  Report  (U)  dated  3Julyl961.  The  complete  report  No.  513 
It  claiiified.  The  purpoieof  this  report  it  to  permit  the  propagation  Infotmatioo  in  Report 
No,  513  to  appear  in  the  unclaitifled  literature. 

The  leductioc  and  analysis  of  the  dau  were  done  by  the  Ctotley  Division  of  Avco 
Corporation  under  Contract  N140(70024)6716SB  with  the  Underwater  Sound  Laboratory. 
R.  W.  Haste,  Jr. ,  and  R.  A.  Wettervelt,  both  of  USL,  provided  technical  cotuulutlons. 
All  the  data  analyzed  were  provided  by  the  Laboratory. 


ABSTRACT 


Project  Colossus  11  wu  established  in  1954  to  Investigate  acoustics  in  shallow  water 
(ISO  lathoms  or  leu).  A  portion  of  that  pro^m  was  devoted  to  a  study  of  underwater 
acoustic  propagation.  The  fretiuetscy  range  of  interest  was  100  to  3000  cps. 

Presented  In  this  report  are  the  results  of  the  acoustic  propagation  loss  studJes.wbich 
include  a  summary  of  the  analysis  and  a  prediction  method  of  estimating  the  stxulstical 
distribution  of  propagation  lou. 


ADMINISTRATIVE  INFORMATION 


Work  on  the  Colossus  II  program  was  accomplished  under  Navy  Project  and  Task 
No.  SS046002-8037  and  USL  Project  No.  1-500-00-00.  This  latter  project  has  been 
terminated  and  the  program  Is  now  Included  under  USL  Project  No.  1-405-00-00. 
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COLOSSUS  n  SHALLOW- WATER  ACOUSTIC  PROPAGATION  STUDIES 


INTRODUCTION 

This  report  covers  the  Colossus  II  study  of  shallow-water  acoustic  propaga- 
tionloss  data.  The  reduction  and  analysis  of  the  data  were  done  by  the  Crosley 
Division  o^  Avco  Corporation  under  Contract  N140{70024)67165B  witli  the  Under¬ 
water  Sound  Laboratory.  R.  W.  Hasse,  Jr.,  and  R.  A.  Westervelt,  both  of 
USL,  provided  technical  consultations.  All  the  data  analyzed  were  provided  ' 
by  the  Laboratory. 

The  bulk  of  the  measurements  was  made  by  the  Laboratory  in  the  period 
since  1955.  Although  coverage  represented  by  the  data  is  incomplete,  the 
entire  collection  now  available  in  the  shallow-water  file  represents  a  degree 
of  completeness  not  likely  to  be  improved  for  many  years.  The  Underwater 
Sound  Laboratory  Colossus  data  are  supplemented  by  data  from  MEDEA, 

ARL,  WHOI,  CUDWR,  and  UCDWR. 

Presented  are  sximmaries  of  the  analysis,  the  status  of  the  shallow-water 
file,  and  a  complete  prediction  method  of  estimating  the  statistical  distribu¬ 
tion  of  propagation  loss. 

No  theoretical  treatment  of  the  data  has  been  made,  the  prediction  method 
being  semi-empirical,  largely  equivalent  to  that  used  in  Pioject  AMOS.  Tha 
validity  of  the  method  depends  on  its  consistency,  reasonableness,  and  utility 
in  representing  the  available  data. 


The  processed  Acoustic  Data  Sources  are  listed  later  in  the  report. 
Other  documents  were  received,  but  owing  to  lack  of  information  they  could 
not  be  processed.  A  complete  bibliography  of  the  Acoustic  Data  Sources 
appears  in  USL  Report  No,  513  (CONFIDENTIAL). 


DESCRIPTIOi;  OF  SHALLOW- 
WATER  PROPAGATION 

A  number  of  formulas  and  charts  are  presented  with  summarize  the 
existing  shallow- water  propagation  data  for  acoustic  frequencies  from 
100  cps  to  2820  cps  and  which  can  be  extended  to  higher  frequencies.  The 
major  source  of  data  for  this  purpose  came  from  the  Colossus  II  pro¬ 
gram  being  carried  out  by  the  Underwater  Sound  Laboratory.  During  the 
years  1956  to  the  present,  large  quantities  of  propagation  loss  data  and 
reverberation  data  were  obtained,  mostly  with  explosive  soiirces,  but 
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also  with  underwater  sirens  and  sinusoidal  acoustic  generators,  off  the 
East  Coast  of  the  United  States.  The  shots  were  recorded  on  magnetic 
tape  and  analyzed  in  logit  frequency  bands  from  70  cps  to  28 ZO  cps. 
Measurements  were  made  at  various  projector  and  receiver  depths  at  ranges 
out  to  50  miles  or  more. 


The  point  of  departure  for  the  present  study  was  the  Project  AMOS  analy¬ 
sis  at  higher  frequencies  in  deep  water  and  recently  completed  studies  of 
surface  scattering  loss  in  isothermal  layers.  The  propagation  loss  data  over 
the  frequency  range  from  100  cps  to  3000  cps  were  analyzed  and 
interpreted  according  to  a  definite  model  compatible  wth  tlie  AMOS 
analysis.  Equations  of  propagation  loss  were  fitted  with  semi -empirical 
coefficients.  The  steps  used  to  arrive  at  this  model  were: 


b.  Studying  acoustic  patterns  in  situations  where  one  of  these  variables 
is  dominant. 

c.  Making  the  simplest  assumptions  regarding  the  acoustic  interrelation 
of  these  variables  and  adding  complications  only  when  necessary  to  incor¬ 
porate  a  large  body  of  data  into  the  model. 

d.  Requiring  that  the  model  be  self-consistent  with  the  AMOS  analysis 
and  the  subsequent  lower -frequency,  deep-water  analysis  based  on  sea- sur¬ 
face  scattering. 


Propagation  Data  Analysis 

The  analysis,  then,  consisted  of  finding  the  most  important  variables  oc¬ 
curring  in  the  shallow- water  propagation-loss  data.  There  was  no  doubt  that 
range  and  frequency  were  the  two  most  important  variables.  The  propagation 
loss  increased  with  frequency  above  100  cps.  One  would  also  expect  the  loss 
to  increase  with  range  following  a  20  log  R  law  out  to  some  range  which 
could  vary  f^om  one  situation  to  another,  and  then  follow  some  different  rela¬ 
tionship  as  boundary  effects  and  absorption  became  more  important.  The 
hypothesis  was  made  to  define  skip  distance  (H)  in  such  a  way  that  H  repre¬ 
sents  the  maximum  range  at  which  rays  first  maice  contact  with  either  the 
sea  surface  or  bottom. 

It  was  shown  in  the  Project  Amos  analysis  thac: 

H  0.5  V  L  upward  refracdoo  in  isotbermal  water  ( 1  A) 

and  _ 

H  0.4  for  downward  refraction  in  aeerafe  ne^tlTC  gradient  water  (IB) 
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where 


L  is  the  layer  depth  in  feet  and 
D  is  the  water  depth  in  feet. 

The  stvidy  was  concerned  for  the  most  part  with  measurements  made 
using  explosives.  For  converting  received  levels  into  propagation  loss  a 
first  approximation  was  made  that  the  best  estimation  of  explosive  source 
levels  might  be  realized  by  assuming  square  law  spreading  to  a  nominal 
range  (usually  3000  yards). 

The  data  were  then  fitted  to  the  form: 

-  20  logjjR  +  *R  +  60  (<n>)  R  <  H  (2A) 

-  10  logjjR  +  10  logjjH  +  aR  +  «'  (R  -H)  +  60  (db)  R  >  H  (2B) 


where 


R  is  the  range  in  kyds, 

a  =  O.Olf*  db/kyd  (average  absorption  coefficiei*'  for  low 

frequencies), 

f  is  the  frequency  in  kc,  and 

aj  is  the  residual  shallow-water  propagation  loss  factor,  db/4cyd. 
(This  varies  Irom  one  situation  to  another.) 

This  varies  from  one  situation  to  another. 


Skip  Distance  as  a  Scaling  Parameter 

The  residual  propagation  loss  factor  a'  is  then  characteristic  of  shal¬ 
low  water  propagation  and  is  a  function  of  frequency,  bottom  type,  sea  state, 
water  depth,  projector  depth,  receiver  depth,  sound  speed  structure,  and 
other  variables.  In  studying  the  depeuuence  of  the  propagation  loss  on  environ¬ 
mental  factors,  it  would  be  extremely  helpful  to  have  a  scale  factor  which 
would  increase  the  effective  number  of  measurements  and  replace  range  and 
water  depth  as  descriptive  variables.  The  skip  distance  (H)  was  found  to 
apply  in  isothermal  layer  propagation  and  was  natural  to  try  in  this  case.  By 
dividing  H  into  the  range  (R)  we  obtain  an  effective  number  of  bottom  and 
surface  contacts  which  control  the  loss  along  a  limiting  ray.  The  propagation 
loss  is  then  controlled  by  those  rays  having  the  smallest  number  of  bottom 
and  surface  contacts.  The  skip  distance  concept  normalizes  the  effects  of 
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the  bottom  and  sea  surface  on  a  **per  contact*'  or  "per  bounce"  basis. 
Thus,  setting  R/H  =  n  in  Eq.  (2B)  we  arrive  at: 


-  201ogjQR  +  aR  +  60  (db) 

■  10  logj^R  +  10  logj^H  +  aR  +  a^  (n  -1)  +  60  (db) 


R  .<  H 

(3A) 

R  >  H 

(3B) 

where 

a^  is  the  shallow- water  attenuation  constant  on  a  per  bounce  basis. 

Actually,  the  transition  from  spherical  to  cylindrical  spreading  should  occur 
over  several  zones.  The  data  indicate  that  this  transition  is  complete  by 
the  ninth  zone.  Taking  an  intermediate  value  for  spreading  during  transition, 
we  have  finally  Eqs.  (4A),  (4B),  and  (4C). 


201og,QR  +  sR  4-  60  (db) 

R  .<  H 

(4A) 

15  1og,,R  +  .P  +  -1^ 

+  5  logj^H  +  60(db) 

8H  >  R  >  H 

(4B) 

101og,,R  +101  +  -1^ 

+  101og,„H  +  64.5  (db) 

R  >  8H 

(4C) 

Effect  of  Thermal  Structure 

The  data  showed  a  strong  increase  in  propagation  loss  in  summer  com¬ 
pared  with  winter.  Since  the  sound  speed  gradients  were  negative  to  the 
bottom  in  the  summer  and  positive  to  the  bottom  in  winter,  we  were  able  to 
compute  quite  definite  skip  distances  for  both  these  situations  from  Eqs. 
(lA)  and  (IB).  For  the  case  where  the  layer  depth  is  less  than  the  water 
depth,  the  following  approximate  expression  was  used  to  compute  the  skip 
distance  of  the  limiting  ray: 


H 


L  D 
8 


kfd . 


(5) 


Eq.  (5)  gives  a  smaller  H  if  L  =  0  than  does  Eq.  (IB)  and  represents 
the  shallow- water  data  better. 
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BOTTOM  ATTENUATION  FACTOR  FOR 
SHALLOW  WATER  PROPAGATION 


Fig,  1  -  Bottom  Attenuation  Factor  for  Shallow  Water  Propagation 


Effect  of  Bottom  Type 

Sines  the  sea  states  were  almost  always  low  in  the  summer,  it  was  possi-  ■ 
ble  to  compare  the  loss  per  bounce  for  sea  state  0  over  mud  and  over  sand. 
These  two  were  the  principal  bottom  types.  Figure  1  shows  the  results. 

It  is  seen  that  the  bottom  loss  for  sand  bottoms  is  appreciably  less  than  that 
for  mud  bottoms,  as  expected.  It  may  also  be  seen  that  these  results  are 
compatible  with  the  deep-water  bottom  losses  as  reported  in  project  AMOS 
for  grazing  between  10®  and  30®. 
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Fig.  2  -  sh*Uow-W*ier  Attenuation  Factor  vs  Sea  SUte  for  Sand  Bottom 


Effect  of  Sea  State 

In  Fig,  2,  the  effect  of  sea  state  is  demonstrated  for  the  winter  data. 
In  this  case  where  upward  refraction  occurs,  it  is  expected  that  the  sea 
state  would  have  a  strong  influence  on  the  propagation.  Actually,  the  loss 
per  contact  looks  like  the  bottom  loss  (most  of  the  measuirements  were 
made  over  sand)  rather  than  like  the  sea  surface  loss  shown  in  Fig.  3. 

The  surface  scattering  loss  has  a  much  larger  frequency  dependence.  For 
example,  this  may  be  seen  by  setting  the  mean  wave  height  equal  to  1  foot 
(sea  state  0-1)  in  the  product  of  frequency  and  mean  wave  height.  This 
behavior  with  sea  state  might  be  expected  if  most  of  the  energy  scattered 
from  the  surface  is  reflected  by  the  bottom  and  returned  for  propagation 
down  the  channel.  In  deep  water,  a,  measures  the  energy  which  is  scat¬ 
tered  ox’t  of  the  surface  channel  and  which  is  either  lost  in  absorption  on  the 
way  to  the  bottom  or  absorbed  in  the  bottom  through  high  grazing  angle 
incidence. 
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Fig.  3  -  Surface  Loss  vs  Wave  Height -Frequency  Parameter 
Surface-Bottom  Coupling 

The  problem  then  is  to  determine  the  extent  of  the  coupling  between  this 
surface  and  the  bottom  from  the  measurements  data  of  Fig.  2  for  sea 
state  0,  sea  states  1  and  2,  and  sea  states  3  and  greater.  If  r,  is  the 
surface  reflection  coefficient,  then  (1  -  r^  )  is  the  surface  scattering  co¬ 
efficient  and  a^  =  -10  logj^  r^  is  the  surface  loss  in  db/bounce.  The  sur¬ 
face  components  and  (1  -  r^ )  must  suffer  different  interactions  with 

the  bottom  if  any  sea  .state  dependence  is  exhibited.  We  would  expect  scat¬ 
tered  rays  to  suffer  greater  loss  because  they  are  steeper.  The  simplest 
expression  satisfying  this  requirement  and  the  observed  data  was: 

^  <1  -  '.)fb 

where 

r,  is  the  fraction  of  energy  tiansmitted  down  the  channel  when  a 
bottom  event  is  coupled  with  each  surface  reflection. 

a,  =  -10  logj^r,  (shallow-water  attenuation  constant,  db/bounce). 
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Tb  is  the  bottom  reflection  coefficient  for  the  angles  of  incidence 
occurring  in  shallow  water  and  is  given  in  Fig.  1  for  sand 
and  for  mud  bottoms.  (3^  =  -10  log,o  r,, ), 

The  observed  behavior  for  a,  has  profound  significance  on  the  inter¬ 
pretation  of  the  shallow-water  propagation  mechanism.  The  fact  that  the 
gracing  angle  rays  represented  by  r^  must  be  multiplied  by  r^  indi¬ 
cates  that  there  is  a  bottom  loss  suffered  by  near-grazing  rays.  The  fact 
that  the  scattered  rays  must  be  multiplied  by  r*  means  that  the  angular 
dependence  of  bottom  reflection  loss  is  such  that  these  steeper  rays  suffer 
twice  the  loss  of  the  near  grazing  rays. 


Interpretation  of  Bottom  Loss  at 
Small  Grazing  Angles 

The  observed  bottom  loss  at  small  grazing  angles  seems  to  imply  a 
mode-changing  process  that  takes  energy  for  a  refracted  ray  system  and 
converts  it  to  a  simple  bottom  and  surface  bounce  mode  at  great  ranges. 
This  is  equivalent  to  the  normal  mode  treatment  of  propagation  in 
shallow  water.  This  mode-changing  process  seems  to  occur  for  any 
sound  speed  structure.  The  sound  speed  structure  determines  the  'kip- 
distaince  zones.  It  would  appear  that  after  the  first  few  zones,  the  skip-dis¬ 
tance  has  a  purely  local  significance  representing  a  limiting  free  path  for 
rays  which  happened  to  be  scattered  into  grazing  angles  at  that  range.  If,  as 
is  shown  in  the  data,  grazing  rays  must  suffer  one  bottom  contact  loss  per 
skip  zone,  then  they  can  not  be  continuous  from  one  zone  to  another.  The 
possibility  that  horizontal  rays  suffer  a  bottom  loss  is  hardly  conceivable. 
The  transition  from  one  propagation  mode  to  the  other  probably  occurs  over 
several  skip  zones^and  we  have  adopted  the  convenient  model  that  it  occurs 
over  the  first  eight  skip  zones.  (See  Eqs.  (4A),  (4B),  and  (4C).  ), 


Near  Field  Considerations 


Results  of  USL  CW  measurements  indicated  that  at  the  short  ranges, 
propagation  was  consistently  superior  to  inverse  square  spreading.  These 
results  are  consistent  with  other  shallow-water  data.  In  addition,  several 
studies  conducted  at  USL  have  shown  average  agreement  between  propagation 
studies  made  with  CW  and  with  explosive  sources.  Thus,  it  would  appear 
that  explosive  source  levels  based  on  published  data  are  satisfactory  for 
converting  from'  received  levels  to  propagation  loss.  The  departure 
of  propagation  loss  in  the  near  field  from  square  law  spreading 
can  be  accounted  for  by  taking  into  account  reflections  in  the 
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Pig.  4  •  Ray  Diagram  for  Direct  Radiation  Zone 


near  field,  using  the  boundary  loss 
values  reported  in  this  study. 

Consider  Fig.  4.  This  diagram 
can  be  expected  to  apply  in  the  region 
reached  by  direct  radiation,  which  is 
essentially  free  from  refraction.  It 
may  be  seen  that  there  is  one  ray 
(1  in  the  figure)  with  no  boundary 
contacts,  two  rays  for  each  order  having  an  equal  number  of  bottom  and 
surface  contacts  (e.g.,  2  and  2');  one  ray  for  each  order  having  one  more 
surface  contact  than  bottom  (3)  and  one  ray  for  each  order  having  one 
more  bottom  contact  than  surface  (3').  The  set  of  rays  is  complete.  Accord¬ 
ingly,  upper  and  lower  limits  on  the  sotmd  field  can  be  calculated  assuming 
incoherent  addition  of  rays  and  the  following: 

1  a  inteasity  of  diiect  rsy 
r  a  surface  reflectioo  coefficieat 


r^  a  bottom  reflectioo  coefficient. 


An  upper  limit  is 


I  +  2l(r,r^  -f  rjr*  +••••) 
+  I(r,  +  r*r^  +.•••) 

+  I  (r^  +  r^r*  +••••) 


IK. 


This  is  an  upper  limit  to  the  field  since  it  neglects  the  excess  spreading 
of  the  higher  orders  of  reflection,  and  changes  of  reflection  coefficient  with 
angle  of  incidence.  A  lower  limit  can  be  obtained  by  considering  only  the 
direct  arrival  and  the  arrivals  having  minimum  angles  of  incidence  at  the 
boundaries.  This  is 

1(1  +  2r,r^  +  ',  + 

The  loss  coefficients  of  this  study  may  be  applied  to  those  results. 

Table  1  shows  the  quantities  =  10  log  and  =10  log  Ky  for  the 
various  conditions  represented  in  the  study. 


9 


Table  1  ' 

NEAR  FIELD  ANOMALY,  1» 


Fiaqiwney 

(cpi) 

1 

Set  Sute  0 

kj,  «  10  logKLt  kjj «  10  log  Ku 

SAND 

Sea  Sute  1  Sea  Sute  2  Sea  Sum  3 

Sea  Sute  4 

Sea  Sute  5 

kL 

ku 

kL 

ku 

kL 

ku 

kL 

ku 

kL 

ku 

kL 

ku 

119 

6.4 

12.3 

6.3 

11.3 

6.3 

11.8 

6.3 

11.8 

6.3 

11.8 

6.3 

11.8 

446 

'  6.1 

10.3 

6.1 

10.3 

6.1 

10.3 

5.9 

9.4 

4.7 

7.7 

4.3 

5.4 

1120 

6.0 

10.0, 

5.8 

9.1 

5.0 

6.7 

4.3 

5.2 

3.9 

4.6 

3.6 

4.1 

2890 

5.8 

9.0  ‘ 

4.1 

4.9 

3.7 

4.3 

3.4 

3.8 

3.0 

3.3 

2.8 

2.8 

112 

6.3 

11.3 

6.3 

11.3 

6.3 

MUD 

11.3  6.3 

11.3 

6.3 

11.3 

6.3 

11.4 

446 

5.8 

9.1 

5.8 

9.1 

5.8 

9.1 

5.6 

8.3 

4.6 

5.9 

4.3 

5.3 

1120 

5.6 

8.2 

5.4 

7.6 

4.5 

5.8 

3.9 

4.6 

3.5 

4.0 

3.4 

3.6 

2820 

4.6 

6.1 

3.7 

4.3 

3.3 

3.8 

3.0 

3.3 

2.7 

2.9 

2.7 

2.9 

There  ie  an  exceedingly  small  amount  of  data  concerning  the  variability 
of  r,  and  ri,  with  angle  of  incidence.  It  is  believed  that  these  reflection 
coefficients  decrease  rather  rapidly  with  increasing  angles  of  incidence.  If 
this  is  true,  the  quantity  IK^  flight  be  assumed  to  be  the  best  representa¬ 
tion  of  the  field  in  the  direct  radiation  zone.  For  comparative  purposes  error 
analyses  are  presented  later  using  both  kj,  and  ky  .  This  near-field 
correction  can  not,  of  course,  be  applied  in  close  proximity  to  the  source, 
since  the  intensity  of  the  reflected  waves  is  overestimated  for  small  separa¬ 
tions  between  source  and  receiver.  For  practical  purposes  the  near  field  in 
shallow  water  can  be  regarded  as  covering  ranges  out  to  about  1000  yards. 
Accordingly,  Eqs.  (4A),  (4B),  and  (4C)  should  be  modified  to  read: 

N,  -  20  log^R  +  sR  +  60  (db)  1  <  R  <  H  (4A)' 

-  15  bg,^R  +  .R  +  +  5  log„H  +  60  -  (db)  8H  >,R  >  H  (4B)' 

-  10  log,„R  +  .R  +  ..,  +  10  log,„H  +  64.5  -  k^  (db)  R  >  8H.  {4C)' 
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OTHER  MODES  OF  PROPAGATION 

The  foregoing  description  applies  to  the  propagation  situation  which  oc¬ 
curs  moat  of  the  time  off  the  East  Coast  of  the  United  States.  Internal  chan¬ 
nels  were  observed  so  seldom  that  they  were  not  considered  an  important 
situation  in  our  analysis.  Presumably,  the  surface  and  bottom  do  not  affect 
this  propagation  mode  which  is  explicable  in  terms  of  spreading  and  temper¬ 
ature  absorption  losses.  This  mode  is  known  to  be  important  in  certain 
localities  such  as  the  Scotian  Shelf.  Very  little  effect  was  found  with  respect 
to  source  and  receiver  depths. 

'No  attempt  was  made  to  carry  the  analysis  below  100  cps.  This  fre¬ 
quency  has  a  wavelength  about  a  quarter  of  the  water  depth  of  200  feet, 
typical  of  these  measurements.  Undoubtedly  a  ray  picture  does  not  apply 
here  and  normal  mode  theory  is  required  to  explain  the  observations. 

Another  phenomenon  associated  with  low  frequencies  which  was  not  treated 
here  was  propagation  from  one  point  to  another  by  way  of  the  sub-bottom  or 
the  seismic  mode.  The  existence  of  this  path  has  definitely  been  establislied, 
and  can  be  of  importance  at  longer  ranges  and  lower  frequencies. 

This  analysis  is  not  applicable  to  bottoms  with  sustained  slopes  in  which 
a  ray  picture  accounting  for  progressively  changing  limiting  angles  must  be 
employed.  The  shallow- water  analysis  applies  up  to  depths  of  100  fathoms, 
which  occur  at  about  100  miles'  from  shore. 

In  the  transition  region  between  shallow  water  and  deep  water,  it  is  known 
that  propagation  may  be  described  quite  well  in  terms  of  a  detailed  ray  analy¬ 
sis.' 


PREDICTION  OF  SHALLOW-  WATER 
ACOUSTIC  PARAMETERS 

A  method  is  presented  below  for  the  estimation  of  the  statistical  distribu 
tions  of  propagation  loss  Nw  in  shallow  water.  Equations  are  given  for 
determining  the  average  values  of  quantities  and  tables  for  the  associated 
problem  errors. 

Mean  Values 

Propagation  loss  is  given  as  functions  of  the  variables  which  follow.  The 
range  of  validity  of  the  variables  is  indicated. 
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(0.1-?.8) 
(100-600) 
(■and  or  mud) 
(0-600) 

(0-5) 


symbols,  the  following  equations  hold  for  mean  value 


N^  •  20  logjjR  +  aR  +  60  (db)  -  k^*  R  <  11 

-  151ogj,R  +  aR  +51og,,H  +60(db)-kj^*  8H  >  R  >  H 

-  lOIog^R  +  aR  +  a^l 

where 

s 

D 


a  r 

+  10  log  +  64.5  (db)  -  k^* 


R  >  8H 


(4A)' 

(4B)' 

(4C)' 


■•VM- 


kyd 


a  >  0.01  f’  approximately. 
•t  ■  “*0‘«8io  f'b'.  + 


or  k. 
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-  1.59  \/A  ,  ft  >  4.14 

-10  log  •  function  of  bottom  type 

-  10  logjg  (1  +  2t,r^  +  f,  +  t^) 


ko 


10  log 
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The  quantities  a,  and  are  displayed  in  Figs.  3  and  5,  re  spec 

tively;  av  is  displayed  in  Figs.  6-11. 
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Fig.  5  -  Bottom  Attenuation  Factor  vi  Frequency 


•TTHIMTIM  »*CTW  lai/Mwaeti 


ATTCMUATIOli  rACTO*  (Aa/MUICCI 


ATTENUATION  FACTOR 
VS  FREQUENCY 

SC*  STATC:  I 
aOTTOM  MATERIALS 

- UUO 

- SAND 


nwouCNCx  t  «cl 


Fig.  7  -  Anenuation  Factor  vs  Frequency  for  Sea  State  1 
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Fig.  8  -  Attenuation  Factor  v$  Frequency  for  Sea  Suie  2 
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Fig.  9  -  Aticnmtion  Factor  vj  Frequency  for  Sea  State  3 
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Fig.  10  -  Anenuation  Factor  vs  Frequency  for  Sea 
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Error  Analysis 


An  error  analysis  was  made  based  on  References  (1),  (2),  and  (10)  in 
Table  4,  p.  27.  since  these  contain  a  large  portion  of  the  data  and  represent  a 
reasonable  cross -section  of  environmental  conditions.  The  analysis  was  car¬ 
ried  out  at  frequencies  of  112,  446,  1120,  and  2820  cps  andrcinges  of  3,  9,  30, 
60,  and  90  kyd.  For  each  document-frequency-range  combination,  approxi¬ 
mately  30  propagation  values  were  picked  at  random.  These  are  called  Nw  . 
i,  Corresponding  to  each  such  value,  avalue  (N^)  was  computed  using  Eqs.  (4A)', 
(4B)',  (4C)',  and  the  near-field  correction  kt,  .  The  differences,  or  anom¬ 
alies  Na  =-Nw  -  Nw,,  were  computed.  Table  2  shows  the  distribution  of 
at  each  range  and  frequency. 


Table  2 

DISTRIBUTION  OF  N^ 


Range  (kyd) 

3 

9 

30 

60 

90 

112  cps 

Upper  Quartlle 

-1 

-1 

5 

7 

C 

Median 

-3 

-4 

-5 

-3 

-3 

Lower  Quartlle 

-6 

-7 

-13 

-10 

-11 

Interquartile  Range 

5 

6 

18 

17 

16 

Number  of  Observations 

67 

72 

74 

70 

60 

446  cps 

Upper  Quartile 

0 

2 

12 

7 

4 

Median 

-3 

-5 

-4 

-4 

-4 

Lower  Quartile 

-7 

-11 

-15 

-14 

-19 

Interquartile  Range 

7 

13 

27 

21 

23 

Number  of  Observations 

92 

107 

105 

107 

93 

1120 

cps 

Upper  Quartile 

2 

5 

17 

9 

5 

Median 

-2 

-1 

-2 

-1 

-7 

Lower  Quartlle 

-6 

-6 

-16 

-12 

-20 

Interquartile  Range 

8 

11 

34 

21 

25 

Number  of  Observations 

92 

108 

103 

108 

78 

2820  cps 

. 

Upper  Quartile 

2 

7 

14 

2 

3 

Median 

-1 

2 

2 

-6 

-9 

Lower  Quartlle 

-6 

-6 

-p 

-17 

-17 

Interquartile  Range 

12 

13 

22 

19 

14 

Number  of  Observations 

90 

100 

97 

86 

70 

Note;  The  above  value!  were  obtained  uilng  the  follovlng  values 
for  kt:  IIJ  7 
446  G 
1120  6 
2820  4 
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Table  3 

DISTRIBUTION  OF  N/ 


Range  (kyd) 

3 

0 

112  cpi 

30 

60 

90 

Upper  Quartlle 

3 

5 

8 

12 

10 

Median 

0 

0 

2 

0 

1 

Lower  Quartlle 

-4 

-3 

-3 

-5 

-6 

Interquartile  Range 

7 

8 

11 

17 

16 

Number  of  Polnu 

45 

84 

446  cpi 

92 

82 

68 

Upper  Quartlle 

2 

4 

7 

13 

6 

Median 

-1 

0 

0 

-3 

-6 

Lower  Quartlle 

-4 

-6 

-10 

-10 

-17 

Interquartile  Range 

6 

10 

17 

23 

23 

Number  of  Poinu 

62 

101 

IICO  cps 

129 

120 

ICl 

Upper  Quartlle 

2 

5 

20 

16 

6 

Median 

-2 

0 

5 

5 

-5 

Lower  Quartlle 

-6 

-5 

-8 

-11 

-17 

Interquartile  Range 

8 

10 

28 

27 

23 

Number  of  Point* 

63 

110 

127 

110 

85 

Reference*  (1),  (2),  (10),  (18),  and  (19) 
Near  Field  Anomaly  »  ku 


A  similar  analysis  was  carried  out  using  the  near-field  correction  factor 
ko  insteal  of  k^,  .  In  this  case  all  available  data  from  References  (1),  (2), 
(10),  (18),  and  (19)  in  Table  4  were  used  for  the  frequencies  112,  446,  and 
1120  cps  and  ranges  3,  9,  30,  60,  and  90  kiloyards.  The  anomalies 
in  this  case  are  defined  as 


where  N*.  is  the  predicted  propagation  loss  using  ku.  The  results  are 
shown  in  Table  6. 

For  112  and  446  cps,  the  results  of  Table  3  show  a  better  average 
fit  of  the  prediction  equations  to  the  data  than  does  Table  2.  For  1120  cps 
the  results  of  the  two  tables  are  comparable.  On  the  average  the  inter¬ 
quartile  ranges  of  Table  3  are  less  than  the  interquartile  ranges  of 
Table  2. 
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Fig.  12  (Coat'd)  •  Distributions  of  Propagation  Lost  Anomaly 

0 

Distributions  of  the  propagation  anomaly,  Na  ,  were  computed  for  the 
frequencies  112,  446,  and  1120  cps  for  the  ranges  3,  9,  30,  60,  and  90 
kiloyards  combined.  These  are  shown  in  Fig.  12.  The  interquartile  ranges 
are  larger  than  might  be  desired.  One  reason  for  this  could  be  the  variability 
in  the  source  level  of  the  explosives  and  another  is  the  inevitable  bias  of 
data  due  to  lack  of  homogeneity.  Thus,  there  are  fewer  measured  values  at 
the  longer  ranges,  because  the  field  is  sometimes  too  weak  to  be  measured. 
Thus,  unweighted  averages  show  a  trend  toward  apparently  unpredictable 
strong  fields  at  the  longer  ranges.  These  fields  are  probably  due  in  part  to 
other  modes  of  propagation,  principaJ^y  seismic.  The  analysis  and  prediction 
data  presented  herein  are  valid  for  values  of  estimated  propagation  loss  less 
than  about  135  db.  Scatter  diagrams  were  prepared  to  show  the  variability 
of  the  measured  loss  with  the  predicted  propagation  loss.  These  appear  as 
Figs.  1 3A,  1 3B,  and  13C. 
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MEASURED  LOSS(N^)  IN  08 


rig.  -  .Scatter  Diagram  —  Predicted  vs  Measured  Loss 


Pig.  14  -  Predicted  vt  Kiutuied  Lon  Uting  Fig.  IS  -  Predicted  vi  Measured  Loss  Using 

for  Slnusolda  1  Sources  for  Sinusoida  1  Sources 


To  further  compare  the  goodness  of  fit  of  the  prediction  equations  when 
using  both  near-field  correction  factors,  Kl  and  ku,  a  scatter  diagram 
was  constructed  of  data  obtained  from  all  available  shallow- water  propaga¬ 
tion  measuremens  made  using  calibrated  sinusoidal  acoustic  generators  as 
acoustic  sources.  These  scatter  plots  are  shown  in  Figs.  14  and  15.  In 
both  the  explosive  results  shown  in  Tables  5  and  6  and  the  sinusoidal 
measurements  shown  in  Figs.  14  and  15,  the  prediction  equations  fit  the 
measurements  better  when  ku  rather  than  ku  is  employed. 
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Fig.  17  -  Summary  Sheet  of  Acoustic  Source  Parameters 
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SviaaM*.  .ry  of  Acouitic  Data 

Th«  acoustic  data  sources  were  available  in  the  form  of  documents  and/or 
IBM  cards.  Table  4  shows  a  complete  list  of  all  sources  received.  Not  all 
of  the  data  contained  in  these  sources  were  usable.  For  this  study,  propaga- 
tion  loss  data  with  a  minimtun  range  of  the  order  of  3000  yards  were  re¬ 
quired,  with  sampling  done  at  1000-yard  increments  thereafter.  The  max¬ 
imum  frequency  of  interest  was  7070  cps.  These  specifications  eliminated 
from  use  three-fourths  of  the  approximately  22,  000  IBM  cards  received. 
For  example,  in  References  (5),  (6),  and  (7)  propagation -loss  data  are 
given  in  range  increments  of  200  to  250  yards.  Hence,  only  every  fifth 
card  was  used.  Some  cards  were  also  omitted  because  the  frequencies  were 
greater  than  7070  cps.  Similarly,  some  cards  from  Reference  (16)  were 
omitted  because  the  minimum  range  was  greater  than  7000  yards. 


Table  4 

UST  OF  PROCESSED  ACOUSTIC  DATA  SOURCES 


Ref.  No. 

No. 

No. 

Average  No, 

Number 

Received  as 

of 

of 

of 

av  uaaavwa 

of  Cards 

Document  (D) 

Rum 

Preq'i 

Ranges/Graph 

or  IBM  Cards 

(1) 

6 

9 

45 

1,501 

D 

(2) 

9 

9 

54 

2,170 

0 

(3) 

2 

10 

60 

239 

D 

(4) 

7 

10 

40 

558 

D 

(5) 

9 

10 

38 

676 

IBM 

(5) 

IBM 

(6) 

8 

10 

45 

673 

IBM 

(7) 

10 

10 

29 

522 

IBM 

(8)  and  (9) 

30 

3 

4 

137 

IBM 

(10) 

10 

8 

53 

2, 178 

D 

(11) 

26 

4 

6 

4,076 

IBM 

(12) 

18 

4 

13 

49() 

D 

(15) 

6 

8 

57 

349 

D 

(16) 

-- 

— 

-- 

-- 

IBM 

(17) 

1 

8 

21 

21 

IBM 

(18) 

6 

6  ■ 

52 

615 

D 

(19) 

5 

6 

47 

862 

D 

Total 

15,067 
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All  the  propagation -loss  data  derived  from  the  acoustic  data  sources 
were  transcribed  into  workable  IBM  Acoustic  Data  Cards,  the  layout  of  which 
is  shown  in  Fig.  16.  Columns  1  and  2  contain  codes  for  identification  of  the 
reference  documents  listed  at  the  end  of  this  section.  Columns  7  to  9  list 
range  in  kiloyards,  while  columns  73  to  76  show  the  common  logarithm  of 
this  range.  Columns  10  to  72  list  propagation  loss  according  to  frequency 
bands.  Columns  78  to  80  provide  identification  of  acoustic  data  in  the 
references.  Colunnns  3,  4,  5,  6,  and  77  are  reserved  for  future  use. 

The  Acoustic  Data  Cards  are  supplemented  with  Summary  Sheets  of 
Acoustic  Source  Parameters.  The  latter  have  been  classified  as  oceano¬ 
graphic,  mechanical,  and  seasonal,  and  are  shown  in  the  Summary  Sheets 
arranged  according  to  runs.  A  copy  of  a  Summary  Sheet  is  shown  In  Fig.  17. 
A  summary  of  the  locations  and  of  the  pertinent  reference  numbers  is  listed 
in  Table  5. 


Table  5 

GEOGRAPHICAL  LOCATION  OF  PROCESSED 
ACOUSTIC  DATA  SOURCES 


Location 

Long  Island  -  Nantucket 
Long  Island  •  Nantucket 
Greenland 

Iceland  -  Faeroe  Islands 
Long  Island  -  Nantucket 
Nova  Scotia 
Fishers  Island 
Block  Island 

South  Carolina  -  Florida 
Western  North  Atlantic 
South  Carolin.i 


Reference 

Numbers 

(1)  and  (2) 

(3)  and  (5) 

Pan  of  (4) 

Part  of  (4) 

(6),  (7),  and  part  of  (10) 
Part  of  (10) 

(8) 

(9) 

(12) 

(15),  (18),  and  (19) 

(17) 
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CONCLUSIONS 


a.  Bottom  loss  in  deep  and  shallow  water  are  compatible.  Differences 
can  be  attributed  to  steeper  angles  of  the  deep  water  rays.  Bottom  loss  over 
mud  is  appreciably  higher  than  over  sand. 

b.  Surface  scattering  is  the  same  for  deep  and  shallow  water. 

c.  There  is  a  strong  surface -bottom  coupling  in  shallow- water  propaga¬ 
tion.  such  tliat  the  propagation  losses  are  controlled  by  the  number  of  con¬ 
tacts  with  both  surfaces.  The  surface  scatters  the  rays,  wheret.s  the  bot¬ 
tom  absorbs  them. 

d.  The  thermal  structure  of  the  water  affects  the  propagation  through 
its  influence  on  skip  distance  and  the  number  of  surface  and  bottom  contacts. 

o.  The  distribution  of  sea  states  in  shallow  water  (100  fathoms)  appears 
to  be  quite  different  from  that  in  deep  water. 

f.  On  the  average,  propagation  loss  is  independent  of  source  and 
receiver  depths. 

RECOMMENDATIONS 

a.  Improve  analytical  model, 

b.  Obtain  a  more  detailed  picture  of  the  loss  mechanisms,  ag  and 

Sa  and  their  interdependence.  Include  analytical  treatment  and  specialized 
measurements.  For  example:  (1)  obtain  a  better  model  of  surface -bottom 
interactions  through  a  transport  type  of  phenomenological  model;  (Z)  make 
direct  measurement  of  bottom  loss  at  small  grazing  angles;  and  (3)  deter¬ 
mine  the  angular  dependence  for  both  as  and  as  . 

c.  To  the  extent  that  further  measurement  programs  are  established,  it 
is  recommended  that  the  coverage  of  environmental  conditions  be  completed 
with  rough  water  measurements  in  summer  and  propagation  measurements 
over  mud  bottoms  in  winter. 
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